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ABSTRACT . 


An investigation was monduckad ieee Che woos Si ba laity 
of visually detecting a dewetted zone in a solid rocket 
Sreeopellant through the use of holographic interferometry. 
Techniques used were simple strain, residual strain and 
time-average holography. As a basis for comparison with 
experimental data a finite element program was used to 


model the test samples. 
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A. DEWETTING 

Dewetting is the phenomenon that makes filled materials 
Meecket propellants) differ from classical viscoelastic 
(unfilled) materials. Propellant consists of a composite 
of solid oxidizer particles imbedded in and bonded to a 
Piel matrix, After an initial linear region, nonlinear 
stress-strain behavior under loading results from the pro- 
gresSive failure of these Pons ae the interface between 
tne fuel matrix and the oxidizer particles. Voids or 
vacuoles are thus created throughout the materials and their 
volume increases with further loading, cauSing an overall 
volume increase in the material. This phenomenon is called 


dewetting. 


ae DEWETTING THEORY 

Lindsey [1] has postulated that a dewetting criterion 
exists and that it 1s expressible as 

J + C4) = GCensetant 

where J5 is the second deviatoric stress invariant, I, the 
Mmifstictress Imvariant and (X¥ 1S a small constant. This 
expression represents a surface in principle stress space 
Poaatets axisymmetric about the hydrostatic axis (C,=0,=03 ). 
Figure 1 is a sketch of the general shape of this surface. 


Poe teomoOreseress. bocatea outside of or on this surface 


Vales cewetted. 





See EVALUATION OF THEORY 

The dewetting created by excess strain is not visible 
to the naked eye. However, changes in the dewetted material 
properties will result in a displacement field that is dif- 
ferent from that of a material which has similar properties 
but does not dewet. To test the postulate of Lindsey, 
holographic interferometry was selected to detect these dis- 
placement changes because (1) it 1s Sensitive to very small 
displacements and (2) holography models need not be trans- 
parent and birefringent as in the case of photoelasticity. 

Attempts to identify a dewetted area through the visual 
interpretation of the interferometric fringe patterns in- 
clude the technigues of small strains, residual strains and 
vibrational time average. 

A finite element program was used as a basis of compar- 
ison between the displacements obtained holographically from 
both dewetting and nondewetting materials of similar shape 


midergoing small strains. 





ith) Uae RIMENTAL APPARATUS 


fee tiae, HOLOGRAPHIC SETUP 

The holographic setup was of a conventional nature. 

An airtable was used consisting of a two-inch thick piece 
Seeplywood, four by eight feet in dimensions, supported by 
four inner tubes resting on a heavy metal table. Very good 
Vibration insulation was obtained. 

A Spectra-Physics 15 milliwatt continuous wave helium- 
neon laser beam was expanded through a spatial filter and 
directed onto the sample. A front surface mirror was 
mounted beside the sample to provide a reference beam. The 
path length from source to sample to holographic plate was 
approximately the same as the path length from source to 
Psermence Mirror to plate. The holographic plate was mount- 
ed directly in front of the sample. Figure 2 is a schematic 
Pmemototograph of the setup. 

The film used was Agfa-Gevaert 8E75 Scientia four-by- 
five inch plates. The developing procedure was five minutes 
in Kodak D-19 developer, 30 seconds in a stop bath and five 


minutes in a fixer followed by a 30-minute water bath. 


Peetu TENSILE TESTING MACHINE 

The use of a conventional tensile testing machine was 
rejected because it has only one moving crosshead. Strain- 
ing a sample from one end results in the center of the 


sample moving about half the distance of the crosshead. 





This movement is undesirable in holographic work where even 
pure rigid body movement causes fringes in double exposure 
interferograms. 

TO minimize the translation of the center of the sample, 
eeecensile testing machine was built that had two moving 
crossheads. Turning the handwheel resulted in the cross- 
heads moving in opposite directions. One complete turn of 
the handwheel resulted in each crosshead moving 0.005 inch. 
mremire 3 1S a photograph of the machine. 

The handwheel had notches in it every five degrees which 
engaged a springwire latch. This allowed the handwheel to 
be turned in multiples of five degrees in the dark for 


adouble exposure work. 


Peet TEST SAMPLES 

The test samples were six-by-six inch squares, one- 
quarter inch thick, with a one-inch diameter hole in the 
center to act as a stress riser. A one-inch piece of wood 
was bonded to two sides of the samples to provide a place 
for mounting in the tensile testing machine. This type of 
geometry was selected for its ease of construction and its 
two-dimensional plane stress properties. 

The sample material was obtained from United Technology 
Corporation, Sunnyvale, California. The unfilled material 
WesmaOcCtee Motor Case liner Of type C-296-182-7 and the 
filled material was rocket propellant of type UTP-3001. 


Figure 4 is a drawing of the sample geometry. 





Pfs oeince MENT FLiELD DETERMINATION 


A. COMPUTER GENERATED DISPLACEMENTS 
1. The Finite Element Program 
As a basis for comparison with the experimental data 
a finite element program called PSELST (Plane Stress Elastic 
Analysis Using Linear Strain Triangles) was used to model the 
test samples. This program was able to handle plane stress 
Problems with up to six different types of material in its 
composition. This is a necessary feature in order to handle 
the material property changes caused by dewetting. The out- 
puts of the program were the stresses and displacements at 
each node in the model. Reference 2 describes PSELST in 
great detail. 
2. The Finite Element Grid 
The finite element grid is shown in Figures 5 and 6. 
Only the upper right-hand quadrant of the sample was modeled 
due to the symmetry of the sample and of the loading. The 
smaller elements were placed around the hole where the 
stress gradients would be the highest. PrewOnagina of the 
meldecGOordinate system was placed at the center of the hole 
with the X-axis as the horizontal lower boundary and the Y- 
axis as the vertical left hand boundary. 
3. The Boundary Conditions 
The boundary conditions on the problem were straight- 
forward. Nodal points on the X-axis were free to slide in 


the plus or minus X-direction only and the nodal points on 





the Y-axis were free to move in the plus or minus Y- 
direction only. Nodal points along the upper boundary 
were all given a uniform positive Y displacement with no 
X displacement to simulate a crosshead moving up and strain- 
ing the sample. The right edge nodal points were a free 
poundary. 
wee tne Untililed Material Solution 

The material properties for the unfilled sample were 
taken to be 470 psi for the elastic modulus and 0.485 for 
Poisson's ratio, which are typical values for this type of 
Material. Experiments have shown that this material is 
linear to about fifteen percent Strain. For this reason it 
was assumed for the computer solution that the modulus of 
elasticity and Poisson's ratio were constants. 

It can be shown that for a displacement boundary 
condition problem, the value of the elastic modulus does 
not affect the resultant displacement field if Poisson's 
ratio is held constant. AS a Sensitivity test, two computer 
Gisplacement solutions were obtained, one with a modulus of 
470 psi and the other with a modulus of 1090 psi. Both had 
peeevetleson S ratio of 0.485. The eepsill eens displacement 
fields agreed to within 0.0003 wave lengths of light for a 
boundary extension of 18.906 wave lengths. 

A second sensitivity test was conducted with the 
elastic modulus held constant at 470 psi and Poisson's ratio 
Peer oo fOr ene £IirSst run and 0.495 for the second. 


The change in the displacement field was small and probably 
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mime the error range of the holographic technique. For 
example, for a boundary extension of 24.136 wave lengths 
of light, the difference between the predicted X displace- 
ments was 0.0579 wave ene one aie Gt Las pOlnt 33% 

The results of the sensitivity tests indicated that 
the computer displacements for the unfilled material were 
unaffected by an elastic modulus change of every element if 
Poisson's ratio was held constant and that the solution was 
relatively insensitive to a Poisson's ratio change of every 
element if the elastic modulus was held constant. 

Peeeebne Filled Material Solution 

The modulus of elasticity for the filled material 
was obtained from a uniaxial tensile test. (Figure 7). 

The stress-strain curve was linearized by dividing it into 
four segments and replacing each segment with a straight 
line. Poisson's ratio for the filled material was estimated 
from data obtained by Lindsey and Wood [1]. Table I isa 
summary of the data. The departure of the stress-strain 
curve from linearity is due to the onset of dewetting and 
the stress at that point, O, Pers Cabiledwtne. yield or de- 
wetting stress. 

The first segment of the curve starting at EGS 
stress was considered the nondewetted portion and was 
designated material of type one. Type two material, which 
was partially dewetted, was represented by the second seg- 
Mememomewthie Curve. Type three was dewetted to a greater 


degree and type four was considered to be completely dewetted. 


uh 





mes upper end Of the fourth segment was the stress, Ores: ; 
at which probable failure of the material would occur. 

To relate the two-dimensional stress of the computer 
solution to the one-dimensional stress-strain curve, an 


equivalent stress o). was used. Equivalent stress [3] is 


defined as 


F 


i, A ie 
J-6 )+(0,-0,) + (6-6) 


as 


O. = SE 
J5= 


on 


Ol 


+6(T cyt T<,+ 1 ) 


BOr plane stress 


Thus 
2. 2 oe 
SG. = VO e020, +37.) 
timeswequivalent stress was used because the value of OG has 
been found to be very small for this type of material. For 
moueal values of CX and the stress field of the specimen, 
the dewetting criterion and the equivalent stress agree to 
within one percent. | 
For the computer solution the upper edge was extend- 
ed 0.01 inches at a time. Although the overall properties 
are nonlinear, each 0.01 inch step was considered small 
enough to be linear. After each extension, the output 
SEDSSses One Oy ana ky at each node were used to determine 


the equivalent stress at that node. If the equivalent 


IZ 





stress was less than O, , the node was considered to still 
be of material type one. When the equivalent stress fell 
within the range of segment two, the node was considered to 
be material of ae two and so on. When the equivalent 
Meee ss excecded Umax, ftailure was probable at that node. 
Each quadrilateral finite element had nine nodes associated 
With it, but only four were nodal points in the grid. The 
program, internal to itself, added a node to the midpoint 
of each side and at the intersection of the element's diag- 
onals. It was decided that whenever any two of the element's 
nine nodes, six in the case of a triangular element, had 
changed to the next higher material type, then the element 
itself would be changed to that next higher material type. 
Once an element or node reached material type four it re- 
mained type four. 

Originally the criterion was that five of nine nodes 
had to change to the next highest material type before an 
element was changed to that type. This resulted pio ar lee ar- hon Rel 
stress buildup at node one and failure was predicted after 
omy 0.2) ianch of upper edge extension. Thus a less strin- 
gent element changing criterion was adopted in an attempt 
to relieve the stress buildup at node one, because it was 
known from experiments that the actual samples did not fail 
Sacre! 0.225 inches of extension. 

Mimo elenmmeercnecking £Or material changes, the 
incremental nodal displacements were added to the nodal 
coordinates to obtain the new coordinates for the next 


extension. 
es 





The incremental procedure used to extend the comput- 
er solution to 0.20 inches may be summarized as the 
following iterative routine: 

eee xeon ENe=upper boundary 0.01 inch, 


b. Add the incremental stresses to a running total 
of stresses from the previous extensions. 


c. Compute the nodal equivalent stresses based on 
the summation of stresses and assign each node 
1ts proper material type based on the linearized 
stress-strain curve. 

dad. Check each element to see if any two of its 
nodes have changed to the next higher material 
type. If yes, change that element's properties 
to the next higher type. 

e. Add the incremental displacements to the nodal 
coordinates to form the new nodal coordinates 
for the next incremental extension. 

feeeketVEn LO Step a. 

Using the above technique, failure was predicted to 
occur at node one after only 0.12 inches of upper edge ex- 
tension. It was known that the samples were extended to 
teezoetnenes before actual failure, so the computer predic- 
ton Of failure was ignored and the solution continued on to 
e220 ianches of extension. 

If the dewetting criterion is applied to the 
@ebascsical solution for the stress field around a hole in 
Mmiiiinite sheet, the outline of the dewetted zone will be 
as shown in Figure 8. For comparison, the computer solution 


for the dewetted zone in the filled material using equivalent 


stress is shown for 0.13 inches of extension in Figure 9, 
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6. Nondimensional Displacements 

The X and Y displacements obtained from PSELST were 
due to a 0.01 inch extension of the upper edge. However, 
the unfilled material was considered linear so that dis- 
placements for extensions other than 0.01 were obtained by 
scaling. 

Predicting displacements for the filled material was 
more restrictive. The displacements had to be scaled from 
the X and Y displacements of the twentieth step of the 
mferemental routine. 

All displacements were nondimensionalized by divid- 


ing by the wave length of light from the Helium-Neon Laser. 


Pee DLSPLACEMENTS BY HOLOGRAPHIC MEANS 
1. Holographic Displacement Theory 

The method presented here was developed by S.K. Dhir 
and J.P. Sikora [4]. If a light-scattering object undergoes 
an elastic deformation between two holographic exposures 
constituting the interferogram, the fringes will be related 
to the displacement field. Let the point P(X,Y,Z) be dis- 
Placed to P' (X+u,Y+v,Ztw) where u, v avi w are the carteSian 
components of the displacement vector (Figure 10). Points 
Pee Yk, 4k) and S({Xom¥o,Zo) are the observatiGn points on 
the hologram plate and the light source respectively. The 
Seraledistance traveled by the light will be (SP+PH) before 


feeemenmteand (SP'+P'H) after movement. The difference in 


i> 





path lengths can be expressed as some multiple of the wave 
Mmemgthi Of the laser light, L, 

nL = (SP'+P'H) - (SP+PH) | (1) 
where n iS any number. If n is a whole number then the 
object point will appear bright because the two waves will 
be in phase. If n is some odd multiple of 1/2 then the 
waves will arrive out of phase and the object point will 


appear dark. 


a 


SP' = [(X-Xotu) 7+ (Y-Yotv) 7+ (Z-Zotw) 2] 
SP = oe revo) 24 (ie20) 21° (2) 
P'H = [(X-Xk+tu) 24 (y-Yk+v) 24 (g-Zk+w) 2)” 
PH = [ (X-Xk) 24 (v-vk) 24 (g-zk) 2]* 
Rearrange equation (1) 
meee’ SP) +(P "H-PH) (3) 


pmiesertuce (2) into (3) and the first term becomes 


oe = so| ee ae 2{(X-Xo)u + (Y-Yo)v + (2-Zo)w] 
(SP) 


h 
jy eye a) 


~~ 


Neglecting small quantities in the expansion of SP' it can 
be shown that 
(ee oP) = _1{ oy 0 ee a) et (Z-Z0) w | (4) 
The eee - PH) can be found by replacing 
e@ensubscript o with k and 1/SP with 1/PH. 


Now 


ny L = Ayu +. BV + CW G5)} 
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Be 





Where 











AL X-XO + X—-Xk 
ney) ee eel 

ey, = Mee) “eed (6) 
SE le al 
Sollee PH 


The subscript k refers to a single point of obser- 
vation on the MoloGeachic plate rgure 11) and ny), 1s the 
respective fringe number on the point being observed. 

The same object point can be observed through three 
different points on the plate. The result is three equations 
meeenree unknowns which can be solved for u, v, and w. 

fitiels mOt POSSIble to determine the location of 
mm@emzero Order fringe, Equation (5) is modified to 

ieee) G = (Ay-Ay 4 )u + Bee ag) mean) me) 
where ety 4 7) 1s the number of fringes that shift through 
Peo ject PpOint going from position k to k+l on the plate. 
When this subtraction 1S carried out, all the terms involving 
1/SP drop out. Thus when using these equations with an un- 
known zero order fringe the location of the source is unimpor- 
tant. This eliminates errors in determining the physical 
coordinates of the source. This advantage is offset by the 
Meeessctty Of combining Ewo observations of the same point to 
determine one equation. The sign assigned to a fringe shift 
BeenoOt important as long as the system used is consistent. 

In theory it should take only three observations to 


determine the three unknowns. However, the value of the 


Ig 





determinant of these equations is very small due to the 
limited size of the plate. With only three equations large 
mmeors Can be expected. The errors can be greatly reduced 
if more than three equations are generated for each object 
point, and a linearly overdetermined least squares method 
is used to solve them. Subroutine LLSQ in the IBM library 
was used. Nine equations were used to determine the three 
unknown displacements. 

The accuracy of the measurements of the physical set- 
Memwes Within 0.1 inch and the accuracy of the fringe count 
[memeccamated tO be + 0.] fringe number. See Table II for 
a sensitivity test. 

A further refinement of Equation (7) is to divide 
through by L. The three ayn then become nondimensional 


displacements. 


Thus, even if the wave length of the laser is unknown, a 
relative displacement field can be obtained. 

The displacement vector obtained from LLSQ may or 
may not have to be multiplied by (-1). The sign of the 
displacement vector cannot be determined from holographic 
measurements alone because the reconstructed image and the 
direction of fringe movement remains unaffected even if the 
direction of the displacement vector is reversed by revers- 
ing the sequence of the two exposures. A priori knowledge 
about the direction of at least one displacement component 
is therefore necessary in order to assign proper signs to 


the other two components. 
18 





If the fringe distortion about a point iS severe, 
Ehere may be larger than 0.1 fringe count errors. A tech- 
nique proposed by Dhir to detect this involves plotting 
Meange difference, (n,-ny4}4) ={\n, along a given cross 
Peetr1on Of the object. A deformed object should show con- 
memuous and smooth fringe differences. If the plots of AG 
are not smooth, a smooth curve is passed through the data. 
in this way the An for a distorted fringe area can be 
obtained with more accuracy. 

2. The Unfilled Sample Experiment 

The sample was extended to 0.135 inches and allowed 
tO relax for more than one-half hour. The first exposure 
was then made and the wheel was turned forty degrees to make 
the second exposure (Figure 12). It was found by experiment 
that for these samples a forty-degree turn of the wheel pro- 
duced the best fringe density. A larger turn resulted in 
mMeewittinges being too close for accurate counting. A smaller 
turn produced fringes that were too broad and increased the 
Gifficulty of estimating the fringe number at a point. 

Nine viewing ports were drawn on the holographic 
Plate with india ink and numbered Snow enrough nine. (See 
Hogure ill). 

The location of the zero order fringe was unknown, 
so two observations of a Single grid point had to be com- 
bined to produce one equation. Each grid point was observed 
nine times to ensure enough over-determination for the 


least squares computer solution. The nine viewing ports 


iL) 





were used in the following combinations for each point: 
m7, 1-8, 1-9, 2-7, 2-8, 2-9, 3-7, 3-8, 3-9 or 1-3, 1~6, 
1-9, 4-3, 4-6, 4-9, 7-3, 7-6, 7-9. If certain combinations 
were unreadable, variations on the above were used. 

For example, when determining the fringe difference 
for set 1-7, the object grid point would be sighted through 
port number one and assigned an arbitrary fringe number. 

If the grid point viewed through port one was half way be- 
tween two dark fringes, it would be assigned a fringe number 
of 0.5. The line of sight was then slowly changed to view- 
port seven while keeping count of how many fringes passed 
over the grid point. If the grid point when viewed through 
port seven was between fringes three and four, then the 
number 3.5 would be assigned. For purposes of the LLSQ 
solution, the fringe difference would be 3.0. The actual 
tenth of a fringe number assigned for each point was deter- 
mined by visual interpolation between centers of dark 
fringes. Most of the fringes could be eee estimated to 
plus or minus one-tenth of a fringe. For crowded or severely 
distorted fringes this was not always possible. This is why 
the fringe-difference curve smoothing was used. 

foeaddttton eo, the 59 grid points used there were 
seven defined points on each crosshead to provide a refer- 
ence for the determination of rigid body displacements and 
the true extension of the sample. 

The smoothed data were used to solve for the total 


i eeoma W Of Gach grid point. The U, V, and W obtained 
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consisted of two parts, the displacement relative to the 


axis system and the displacement due to rigid body affects. 


Tus : 
U = Ud + Urb 
V = Vd + Vrb 
W = Wd + Wrb 


The subscript d refers to displacements relative to 
Pewaxes. This 1S the quantity desired for comparison with 
Ppewcomputer solution. The subscript rb refers to the rigid 
body displacement. The rigid body displacements were deter- 
menecaeand Subtracted from the total displacements to yield 
the relative displacements. 

The rigid body displacements of the sample were 
determined from the displacement of the X and Y axes. It 
was assumed that the X and Y axes remained straight even 
under stress due to the geometric symmetry. Therefore, any 
Peer oplacement along the X-axis was due completely to rigid 
body displacement. Any U displacement along the Y-axis was 
due completely to rigid body displacement. That is, U = Urb 
along the Y-axis and V ne ee rong tie) X-axis. 

As stated previously, a priori knowledge of one 
Meelacement component must be known before the proper sign 
Can be given to the other two components. It was known that 
the V displacements of the nodal points on the Y-axis were 
positive since that was the sign of the upper crosshead 
extension. Therefore the sign of the U and W displacements 
was also determined for the Y-axis. A linear least Squares 


curve was fitted through the U displacements, making Urb a 


Ze 





meumettOon Only of Y. It was then assumed that the Urb dis- 
Sacement Of any point on the body was a function of its Y 
coordinate only. | 

By letting Y equal zero, the Urb for points on the 
X-axis was determined. Although Ud along the X-axis was 
unknown its sign was known to be negative from the Poisson 
effect of a positively strained body. Therefore, 

ie Urb = 0 
and in addition Ud must become more negative as the X 
distance from the origin increases. 

As an example, the following is a sample of the data 


Serained £LOr two nodal points on the X-axis. (Figure 6) 


Point U We W 
3 Bee 52 S57 45 SO o2) 
6 Ldn 2S Paes 3 Oh PSV OS, 


All the values are positive but it was not known whether 
they should be. It was determined using the method described 
meres Urb £Or the X-axis was ~2.5057. Applying the formula 
to point three, | 

ome Si eae= 2.5057) = 0 
Shows that the U for point three must be negative and there- 
fore U and W are also negative. Then Ud3 = ~-5.7455 and 
Ud6é = ~8.7658 and the second requirement 1s met. With the 
Vrb for the X-axis known to be negative, these points were 
also fitted with a linear least squares curve and it was 
assumed that the Vrb of any point was a function of its X 


Seeordinate only. 
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iMiemeeterininatton Of Wrb 1S very difficult and can 
only be approximated at best. There are no sections of 
the sample that can be assumed to remain straight for dis- 
placements in the Z direction. In addition, the maximum 
Wd displacement was predicted from the computer solution to 
be only on the order of two wave lengths with most points 
being on the order of one or less, which borders on the 
expected accuracy limits of the fringe count displacement 
prediction method. For this reason W was not considered 
in the comparisons of displacements. 

AS a point of interest, it was noted during the 
counting of fringes that the maximum fringe shift would 
occur if the observation viewports used were ina line 
parallel to the maximum displacement component. 

For the unfilled material experiment the additional 
extension of the upper end was found to be 24.7041 wave 
lengths of light. Forty degrees of handwheel turn was 
expected to cause 22.3 wave lengths of additional extension. 

3. The Filled Sample Experiment 

The methods used to make the interferograms and 
obtain the data for the filled sample were identical to those 
of the unfilled sample. The only change was in the amount of 
extension given the filled sample prior to the making of the 
interferogram. It was found by experiment that the samples 
would fail by cracking when the extension was about 0.225 


inches. The sample was extended to 0.20 inches to ensure 
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that dewetting would take place. Figure 13 is a photograph 
of the fringe pattern obtained. 

For the filled material experiment the additional 
extension of the upper edge was found to be 18.9060 wave 


lengths of light. 
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Vee oLiei METHODS 


A. SaMPLE EXTENSION 

The simplest method attempted was to apply a small ex- 
tension between exposures and visually search for anomalies 
in the fringe pattern. No such pattern was detected which 
could be attributed to a dewetted zone. Figure 13 is a 
typical example of a fringe pattern obtained. Unwanted 
rigid body motion may have caused fringe patterns which 


disguised the zone. 


Peme RESTDUAL STRAIN 

The residual strain method makes use of material prop- 
erty changes. The object was photographed at one state of 
strain and then the crosshead was extended another 0.005 
inch. The increased strain caused more of the sample to 
dewet. The crosshead was then returned to its original 
position and the second exposure made. The nacaniadl ehat 
changed properties between exposures will not be able to 
return to its original state and this will cause fringes. 
Figures 14 and 15 are examples of this technique. The 
resulting fringe patterns vary widely. On several occa- 
Sions, patterns appeared in one or more quadrants that may 
Or may not have been dewetted zones. Fringe patterns on the 
crossheads indicated that they were not returned to exactly 
mg@enwOrtbgGinal extension and/or that rigid body motion had 


taken place between exposures. 
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eee LiME AVERAGE 

Samples of the filled material were acoustically excited 
and time average holograms were made. The range of ex- 
Citation frequencies used was from 60 to 1000 cycles per 
second. The holograms obtained showed no indication of a 


G@ewetted area. (Figure 16) 
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V. RESULTS AND CONCLUSIONS 


Pee tue UNFILLED MATERIAL 

The computer-predicted displacements and the measured 
displacements for the unfilled material are plotted versus 
radial distance for each ten-degree radial of the finite 
element grid. (Paagces 17 through 21). An examination of 
the plots indicates that the linear theory does not accu- 
rately predict the behavior of the sample. The measured 
U and V displacements are higher than predicted on every 
maaial, 

There are indications that the deviation from linearity 
depends upon the intensity of the stress. It can be seen 
that the V measured points match very well with the shape 
of the V predicted curve at the higher radials. Since the 
displacement at any point is the integral of the strains 
femeo that point, it would indicate that the material at 
the higher radials 1s behaving linearly and the larger than 
predicted displacements are due to the nonlinear strains on 
the lower radials. 

Another unexpected result is that the measured U for 
the ten-degree radial is larger than that of the zero-degree 
Paatal. The U for the 20-degree radial is smaller than the 
Ueeor the 10-degree aaaied Pitesti llarger than the zero. 
It is not until the 30-degree radial that the measured U is 


smaller than that of the zero-degree radial. The predicted 
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U is greatest for the zero-degree radial and decreases for 


Seen increase in radial. 


fee los FILLED MATERIAL 

BOolethe filled material the linear predictions, the 
dewetted predictions and the measured displacements are 
plotted versus radial distance for each radial. (Figures 
meeenrough 26). 

The linear V displacements predictions were in general 
better than the dewetted predictions. This 1s due in part 
to the early prediction of failure from the dewetted pro- 
gram. The failure prediction was discounted and the result 
was that the predicted dewetted zone grew to cover most of 
the right-hand side of the grid. The lowered elastic moduli 
of the dewetted lower radials allowed too much strain to 
occur in the elements along the X-axis, with the resultant 
larger displacements at the higher radials. 

Both computer methods underestimate the U displacements 
and they agree with each other very well for the 30-degree 
radial and miloneks. 

The measured values of V fell between the linear and 
dewetted computer solutions but were closer to the linear 
Seeieion. This indicates that the material did dewet but to 


a much smaller extent than predicted by theory. 


ae tbe DISPLACEMENTS COMPARED WITH UNFILLED DISPLACEMENTS 
The measured displacements for the filled and unfilled 


material are plotted versus radial distance for each ten 
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Smeeee racial (Figures 2/7 through 31). The unfilled dis- 
Mereements were converted to the same upper edge extension 
as the filled material by multiplying by 


oO = 


The unfilled displacements are greater than the filled 
displacements for every radial but in general the shapes 
of the displacement curves are Similar. This appears to be 
due to nonlinear straining of the unfilled material in the 


ibeorwer radials. 


ieee ths HOLOGRAPHIC DISPLACEMENT DETERMINATION METHOD 

The method for obtaining general displacements from a 
holographic interferogram proved to be simple to understand 
and easy to use. However, its use is not recommended for 
general use because of the large amount of time that must 
Bemsopenc to obtain data. It 1S estimated that a total of 
twelve hours were spent just to determine the nine fringe 


Gifferences for each of the 73 points used in the analysis. 


E. THE TENSILE MACHINE 

emcee machine, as designed, was not rigid enough. 
Wore! unwanted io CtoMemOuschoemacnine Inducca fringe patterns 
onto the sample that disguised the true patterns and made 
Simple viSual inspection of the patterns impossible. Any 
continuation of this work should be done with a much more 


meera tensile testing machine. 
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fe BOULVALENT STRESS 

The concept of equivalent stress was derived for use 
in predicting plasticity in metals. It is used to relate 
a three-dimensional stress state to a one-dimensional one 
such as a stresS-Strain diagram, however it is only a 
Mmeeemlate, not a proven fact. The effects of plasticity 
and dewetting, as far as the stress-strain diagrams are 
concerned, are the same so it seemed a logical extension 
of the equivalent stress idea to apply it to dewetting of 
meopellants. The results indicate that this is not the 
case and that equivalent stress may not be applied to 


dewetting with accurate results. 


fee lHe FINITE ELEMENT PROGRAM 

From the inaccuracy of the computer displacement pre- 
dictions and the early prediction of propellant sample 
failure it is evident that a finite element material 
Pee orey Change Criterion based on equivalent stress alone 
will not Givemwecaeilshactony results. The effects of visco- 
elastic stress relief will have to be incorporated into the 
computer dewetting criterion so that stresses will build 
more slowly without over-predicting the resultant 
displacements. 

The linear finite element program more closely predicted 
the V displacements of the propellant but it also predicted 


a stress of over 200 psi at node one for a 0.20 inch extension, 
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which is over twice the failure stress obtained from the 
tensile test. The sample did not fail at 0.20 inches of 
extension which emphasizes the importance of viscoelastic 


stress relief incorporation into the finite element model. 


H. CONCLUDING REMARKS 

The comparison of the computer predicted displacements 
to the measured data was unsuccessful due to the inconsist- 
ency of the results. Modifications to the computer tech- 
nique as discussed in section G should produce better 
results. In addition, a more rigid tensile testing machine 
would eliminate errors resulting from rigid body motion and 
unsymmetric crosshead extension. 

The dewetted zone was not detected visually. The most 
promising technique appears to be the residual strain method 


and further investigation of the technique is recommended. 
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LAG 1 


TABULATION OF MATERIAL PROPERTIES 


MATERIAL TYPE POISSON S RATIO ELASTIC MODULUS 
ih 0.485 HO 0 
Zz 0.483 205i 
3 Oa 1 (eo 
4 0.480 nO Aves 
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AIEEE peo 1s 


feof OoL SENSITIVITY Of LEAST SQUARES FRINGE DIFFERENCE 


Data used: 


METHOD 











Object Number First View Port Second Port (Nye 
1 1 3 Se 7 
Ik aL 6 Sec 
il a 9 Cue l 
1 4 2 4.0 
1 4 6 Avil 
1 4 9 Se 
1 7 3 Ae 
1 y 6 4.4 
1 7 9 Bre 
Test results: 
Number of Solution 
Equations used W) V W 
4 i635 -3.172 -19.69 
5 at Ano oo -21.83 
6 -7.875 -~3.020 -17.27 
v7 -~7.894 -2.489 ~14.38 
8 -7.896 -2.542 lar Oy 
9 -7.829 -~2.544 S210 7 
9 -7.799 -2.550 oom Oce mi) 
9 -8.053 -2.544 -15.08 note (2) 
9 -7.605 ~2.544 -15.06 note (3) 


Note 
iecOm 


Note 


Note 


ow 


Odd numbered data cards had 0.1 subtracted 


n and even numbered cards had 0.1 added. 


i) 
3) 


Acco meicemhad 0. i) added to /\n. 


All nine data cards had 0.1 subtracted from ine 
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General Shape of Dewetting Criterion Surface in Principle 
stress Space 
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Figure 2 Holographic Setup 
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Tensile Testing Machine 


Figure 3 
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Test Sample Dimensions (in) 


Figure 4 
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Finite Element Grid with Nodes Numbered 


Figure 6 
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Figure 8 
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Finite Element Grid Showing Predicted Dewetted Zone in 
Filled Material at 0.13 Inches Extension 


Figure 9 
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Schematic for Theoretical Analysis 


Figure 10 
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Holographic Plate with Nine Viewports Drawn on the Glass 
Wren indaa nik 


Reconstructed 
Neledrtaymuc Image 
Plate 





Method of Fringe Counting 
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Time Average Hologram of Filled Sample Excited at 120 CPS 


Figure 16 
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Unfilled Nondimensional Displacements vs Radial Distance 


Figure 18 
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Unfilled Nondimensional Displacements vs Radial Distance 


Figure 19 
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Unfilled Nondimensional Displacements vs Radial Distance 


Figure 20 
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Unfilled Nondimensional Displacements vs Radial Distance 
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Filled Nondimensional Displacements vs Radial Distance 
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Filled Nondimensional Displacements vs Radial Distance 


Figure 24 
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Filled Nondimensional Displacements vs Radial Distance 
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Filled Nondimensional Displacements vs Radial Distance 


Figure 26 
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Nondimensional Displacements vs Radial Distance 
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Nondimensional Displacements vs Radial Distance 


Figure 28 
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Nondimensional Displacements vs Radial Distance 


Figure 29 
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Nondimensional Displacements vs Radial Distance 


Figure 30 


62 








Git, alae! 
6 U unfilled 
ery £iiled 
©) vV unfilled 





-D 


Nondimensional Displacements vs Radial Distance 


Figure 31 
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